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Abstract IC 4406 is a large (about 100'' x 30'') south- 
ern bipolar planetary nebula, composed by two elon- 
gated lobes, extending from a bright central region, 
where there is evidence for the presence of a large torus 
of gas and dust. In this poster we show new obser- 
vations of this source performed with IRAC (Spitzer 
Space Telescope) and the Australia Telescope Compact 
Array. Although the possibility for faint extended emis- 
sion to be missing in the radio maps cannot be ruled 
out, flux from the ionized gas appears to be concen- 
trated in the bright central region. Comparing ATCA 
to IRAC images, it seems that, like in other planetary 
nebulae, ionized and neutral components spatially co- 
exist in IC 4406. 
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^ 1 Introduction 

IC4406 is a well studied southern planetary nebula. It 
has been imaged with several telescopes and at different 
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wavelength ranges. It shows two H2 lobes [6], orthogo- 
nal to the nebula's major axis and each 15" away from 
the center. These peaks are approximately coincident 
with the two blobs observed in Ha+NII and OIII [4], 
interpreted as indicative of the presence of a dense equa- 
torial torus of dust. CO maps [4J show the presence of a 
collimated high velocity outflow in the polar direction. 
Hubble images in Nil, Ha and OIII have revealed the 
existence of an intricate system of dark lane features 
[1], which led to the name of "Retina Nebula" for this 
object p. 

We have observed this source in the radio range to in- 
spect the distribution of the ionized gas in its envelope 
and in the infrared to check for emission from the equa- 
torial dust and molecular gas. 

2 Observations and Data reduction 

Radio observations were performed at the Australia 
Telescope Compact Array on November 24, 2005 (UT: 
from 17:00:00 to 08:00:00) and December 11, 2005 (UT: 
from 15:30:00 to 02:00:00), simultaneously at 4.8 and 
8.6 GHz. The November run was performed with the 
array in 1.5C configuration, while for the December one 
the configuration was 6. OA. The adopted configurations 
are both linear but with different antenna positions, 
giving maximum baselines of 4500 m (1.5C) and 5939 
m (6. OA), minimum baselines of 77 m (1.5C) and 337 m 
(6. OA). The pre-calibration of the array was performed 
observing 0823-500, while the absolute flux calibrator 
was 1934-638. Another target was also observed during 
our two runs and the total on-target time was about 7 
hours for each of the two. The phase calibrator chosen 
for IC 4406 was 1431-48, which is 4.76° away from the 
target. The data were reduced with MIRIAD, follow- 
ing a standard reduction procedure as recommended in 
the MIRIAD User's Guide. The data from the two runs 
were combined into one dataset, obtaining a t^i' coverage 
from 0.9 to 96 kA at 4.8 GHz and from 1.5 to 172 kA at 
8.6 GHz. These correspond to an angular resolution of 
2.2" at 4.8 GHz and 1.2" at 8.6 GHz, while the largest 
observable structures (Largest Angular Scale) are 230" 
and 140" respectively. 

Infrared observations were performed with the Infra- 
Red Array Camera onboard Spitzer Space Telescope 
at 3.6, 4.5, 5.8 and 8.0 /im on March 06, 2004 (UT 
09:54:16.311). Six High Dynamic Range 30 sec dithered 
frames were obtained at each wavelength, for a total ex- 
posure time of 180 sec per channel. In TablHthe point- 
ing coordinates for the two datasets are listed. 

Basic Calibration Data were retrieved from Spitzer 
archive, cleaned to correct such artifacts as mux-bleeding 
and banding and then coadded using IRACproc [5j. 



2 



Table 1 J2000 pointing coordinates for our target and its phase 
calibrator, used for our ATCA run, are shown. 





RA 


DEC 


SST 






IC 4406 


14:22:27.66 


-44:09:00.0 


ATCA 






IC 4406 


14:22:26.28 


-44:09:00.0 


1431-48 


14:35:16.80 


-48:21:47.76 



Table 2 In hues one and two the wavelengths and relative flux 
measurements from Spitzer Space Telescope are listed. Lines 3-5 
summarize the radio observations: frequencies and relative flux 
values in columns 3 and 4, the rms in column 5 (the latter is the 
same within errors for both frequencies). Flux unit is mJy. 



SST 


3.6 jim 


4.-5 jim 


5.8 jim 


8.0 jim 


IC 4406 


62.7 


110.0 


132.4 


307.2 


ATCA 




4.8 GHz 


8.6 GHz 


(7 


IC 4406 




96.47 


91.15 


0.04 


1431-48 




1020 


650 


1 



The whole nebula in each image was boxed with a poly- 
gon and the flux within such region summed up. The 
same procedure was adopted for the field star observed 
West of the central core, so that its fiux was subtracted 
to the overall nebula's one. The result was corrected for 
an infinite aperture, according to the IRAC Data Hand- 
book and then converted into mJy (SST BCD files are 
in units of MJy/sr). 

Tab 12] summarizes the results of our radio and infrared 
observations. 



3 Results and Discussion 

Our radio maps (FigH]) show the presence at 4.8 GHz 
of a 42" X 56^' (3cr level) emitting region, elongated in 
E-W direction; at a 10% of the peak level the size of the 
emitting region is restricted to 32'' x 32''. At 8.6 GHz 
the 10% of the peak level gives a size of 36" x 40", in 
agreement with its 3<j level size. The maps don't show 
any N-S blobs of emission. What is shown is a very intri- 
cate system of emitting lanes (Figl2j), which resembles 
what is observed by Hubble Space Telescope [Ij. An 
inspection of the amplitude vs. nt^distance plots leads 
to the conclusion that faint extended emission may be 
missing in our radio maps despite the combination of 
an extended and a compact array configuration. Figl3] 
shows IRAC images of the central equatorial area of 
the nebula. They were plotted with a linear scale hav- 
ing the peak fiux and 50% of it as thresholds. Channel 

4 is clearly showing the emission from the torus of dust 
surrounding the central star. Its size is about 28" x 20", 
elongated in N-S direction and the angular distance be- 
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Fig. 1 Radio maps of IC 4406 {top: 4.8, bottom: 8.6 GHz). 
The shown contours indicate the zb3cr levels. Both maps were 
produced with NATURAL weights. The synthetic beam is shoiwn 
in the bottom right of each map and the flux density unit is 
Jy/beam. 




Fig. 2 Zoom of the 8.6 GHz map. This map was obtained with 
ROBUST weighting, to balance better resolution and sensitiv- 
ity to faint structures. The RESTORed map was scaled to 1.5" 
beam, shown in the bottom right. The flux density is in Jy/beam 
and the plot ranges from the 3a level to the peak value. 
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Fig. 4 Images obtained with IRAC. From top to bottom: overlay 
of IRAC channels (blue is 3.6, green 4.5, orange 5.8 and red is 8.0 
/im); as in the previous image but not including 5.8 /im channel, 
for better viewing the central core (torus) emission; 5.8 /im image 
scaled for a better view of the arch structures mentioned in the 
text, some of which are pointed by the arrows. 

Table 3 Parameters adopted to match observational points in 
Fig|5] Emission measure is in cm~^- pc, temperatures in Kelvin. 



IRAS IRAC 2MASS Central star 
T=72 T=490 T=1400 T=10^ 



Radio 

EM=9.696 • 10^ T=10000 
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Fig. 3 Zooms of the central core emission observed with IRAC. 
The scale is linear: maximum is the emission peak in the area, 
minimum is 50% of the maximum. Clockwise from top left: 3.6, 
4.5, 5.8, 8.0 fim. 



tween its peaks is about 14^'. The overall size of the 
torus matches with the approximate size of the neb- 
ula in N-S (30''), although its peaks are much closer 
to the center than the H2 blobs repoted in [6], which 
can indicate the torus is partly shielding the molecu- 
lar gas from UV radiation. FigH] is a combination of 
IRAC channels in linear scale. Despite the lower reso- 
lution compared to Hubble images, IRAC is able to de- 
tect the faint emission from the neutral components in 
the envelope and reveals the structure of the elongated 
lobes. They show arches at different distances and in- 
clinations from the central star, connected to the mass 
loss history. The arches which are angularly closer to 
the central star show larger blue emission. This implies 
they have higher temperature, which can be explained 
if they are intrinsically closer to the star. The overall 
structure observed in the envelope corresponds to the 
idea of the central torus as a main collimating agent. 

Collecting literature data and adding our new mea- 
surements, we have built up the SED for our target from 
the radio to the near-IR. We have overplotted to the 
data the sum of the theoretical curves for free- free emis- 
sion and blackbody radiation. Three blackbodies with 
different temperatures were used to match the infrared 
data. The radio free-free emission was reproduced ac- 
cording to the equations in ^. One more blackbody 
was plotted to match one optical point, to correct near- 
IR emission at shorter wavelengths. TabjS] summarizes 
the parameters used for the plot. The curves were sep- 
arately normalized to one point before being summed 
up to match the data points (Fig (5]). As reported in [4], 
more than one dust component is necessary to match 
the observations. Since fewer points were used in that 
study, a different number of components was found. 
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Table 4 Parameters calculated for two reported distances to our 
target. Density and ionized mass are calculated from 5 GHz data, 
according to [2 , dust mass from 25 /im data, according to [3 . 



Distance (kpc) 


0.7 


1.6 


'^ionized (Mq) 


0.11 


0.88 


Mdust (10-5 Mq) 


4.54 


23.7 




4.06 • 10-4 


2.68 • 10-4 


Density (cm~^) 


1.37 • 103 


0.90 • 10^ 



IC 4406 




Fig. 5 Spectral energy distribution with literature (solid circles) 
and new (open squares, with errors within each square) data. The 
single contributions are also plotted. Long dashed curve: free-free; 
dashed: blackbody at 72 K; dash-dot: blackbody at 490 K; dash- 
dot-dot: blackbody at 1400 K; dots: blackbody at 10^ K. 



Although the possibihty of H2 hne emission in IRAC 
bands cannot be ruled out, thus changing the contin- 
uum shape, according to our calculations, cold, warm 
and hot dust need to exist in the circumstellar envi- 



ronment to explain current observations. Tab HI reports 
some parameters calculated assuming two different dis- 
tances as reported in the literature ([^, 

We can speculate that in such a diversified dust en- 
vironment further lower temperature components may 
exist and high sensitivity, high angular resolution obser- 
vations with ALMA will sure give a fundamental con- 
tribution to understand the physics of circumstellar en- 
velopes in planetary nebulae. 
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